function [ ] = plotStrain( domain2D, rx, ry, strain )
%PLOTRESULT Summary of this function goes here
%   Detailed explanation goes here

dim = size(rx, 1); 

px = [domain2D.x'; rx];
py = [domain2D.y'; ry];

n = size(px, 1);

c = zeros(n, 1);

for i = 1 : dim
   vc = max(eig(strain(:,:,i)));
   c(i+n-dim) = vc;
end

c

figure;

tri = delaunay(px, py);

trisurf(tri, px, py, c, c, 'facecolor', 'interp', 'edgecolor', 'r');

colormap(jet(256))

view(2);

axis ( [ -1, 1, -1, 1] );
% caxis( [ 0, 1.7] );
axis square;

colorbar;

end

